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Abstract 

Recently, wireless communication industries have begun to extend their services to machine-type 
communication devices as well as to user equipments. Such machine-type communication devices as 
meters and sensors need intermittent uplink resources to report measured or sensed data to their serving 
data collector. It is however hard to dedicate limited uplink resources to each of them. Thus, efficient 
service of a tremendous number of devices with low activities may consider simple random access 
as a solution. The data collectors receiving the measured data from many sensors simultaneously can 
successfully decode only signals with signal-to-interference-plus-noise-ratio (SINR) above a certain 
value. The main design issues for this environment become how many data collectors are needed, how 
much power sensor nodes transmit with, and how wireless channels affect the performance. This paper 
provides answers to those questions through a stochastic analysis based on a spatial point process and 
on simulations. 

Index Terms 

M2M, stochastic geometry, spatial reuse, outage probability, network design, Poisson point process. 



I. INTRODUCTION 

Wireless personal communication enables ubiquitous exchange of various data types such as 
voice, video, photos, and text among individuals. The emergence of new advanced systems such 
as the IEEE 802.1 lac HI and the 3GPP LTE- Advanced are expected to achieve additional 
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data rates. Of late, wireless communication industries have begun to discuss their scenarios 
serving machine-type communication devices such as meters/sensors as well as user equipments 
such as smart phones BUflH. These machine-to-machine (M2M) communications have extensive 
applications, from monitoring environments to full electrical/mechanical automation (e.g. smart 
grid, smart city, Internet of things), which has been being considered as one of the most crucial 
technologies in future [|5]|[|6]|. The sensor network can also be regarded as a kind of M2M, and 
there have been many studies in the form of ad-hoc networks [JZljll. This paper only considers the 
environment with specific data collectors directly communicating with sensors. This environment 
is suitable when sensor nodes support only simple single-hop communication functionalities and 
deployment of many data collectors is easy. This type of M2M communication is similar to cel- 
lular communication systems where the base stations serve user equipment within their coverage, 
but it has the unique characteristics [ED[|9]: there can be a huge number of devices (e.g. trillions) 
each of which has only a small amount of data and a low activity, and their functionalities have 
to be simple. These characteristics may require technologies differentiated from the conventional 
high data rate human-to-human (H2H) communications. For example, machine-type devices such 
as meters and sensors need uplink resources intermittently for reporting measured or sensed data 
to their serving data collector, but it is hard to dedicate limited uplink resources to each. Thus, 
simple random access can be considered as a solution for directly transmitting measured data 
or initially requesting uplink resources. The data collectors that receive many sensors' measured 
data simultaneously can successfully decode only signals with signal-to-interference-plus-noise 
ratio (SINR) above a certain value. In order to keep a high success probability of many sensor 
nodes' intermittent transmissions, the system may need a lot of data collectors, and conventional 
macro/micro base stations may not be appropriate for these roles. In other words, data collectors 
have to be easy to deploy and cost-effective. They support only simple functionalities and are 
interconnected with external networks through wired or wireless links. It can be considered 
that not only a new type of device for data collection is defined but also such devices as 
pico/femto base stations around sensor nodes play the role of data collectors. Fig. \T\ shows a 
system architecture with data collectors and sensor nodes. In this environment, some questions 
are: How many data collectors are needed? How much transmit power sensors have to use for 
successful transmission? And, how the wireless channels affect the performance. This paper will 
provide answers to those questions through a stochastic analysis based on a spatial point process 
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and on simulations. 

The main factor of determining system performance is the interference from neighbor sensor 
nodes. This interference depends on the spatial distribution and sensor-node access methods. 
Because the spatial configurations of transmitting and receiving nodes can have enormous pos- 
sibilities, it is impossible to consider each possibility. Stochastic geometry provides a useful 
mathematical tool to model network topology, and it also enables analysis of essential quantities 
such as interference distribution and outage [flQl [ f TTl[fl2l. This stochastic geometry has mainly 
been applied to pure ad hoc networks and their performance has been analyzed under the 
assumption of random transmitter location and receiver with fixed distances to its transmitter 
ifTOl lfT3l |f l4 j . This paper considers the environment where both transmitters (sensor nodes) and 
receivers (data collectors) are randomly deployed and each transmitter are served by the data 
collector nearest to it. lfT5l - lH9'l have analyzed the distribution of signal-to-interference ratio 
(SIR) or SINR in random cellular networks where both transmitter and receiver are randomly 
located; 0/7] analyzed the distribution of SIR considering the path loss and shadowing, lfT8l 
derived a simple-form SINR distribution in case of Rayleigh fading and a path-loss exponent of 
four, and [|T9l expanded the analysis results in [[181 into the results for a more general fading 
model including Nakagami-m fading. But, they assumed that each base station always has the 
user equipment within its coverage and communicates with a user equipment that is scheduled 
exclusively within one cell and focused on transmitter-centric coverage (i.e. downlink). [|20l 
modeled CDMA uplink interference power as a log-normal distribution using the moment- 
matching method. Also, [12711 asymptotically analyzed uplink spectral efficiency in spatially 
distributed wireless networks, where the base stations have multiple antennas, by using infinite- 
random-matrix theory and stochastic geometry. The current system is similar to the uplink cellular 
systems, but this paper will only consider random access without any explicit scheduling for the 
simple functionalities of sensor nodes and data collectors. 

The three contributions of this paper are: First, an analysis shows how the channel affects the 
SIR distribution for Nakagami-m fading. A simple form on the SINR distribution is found for 
some special channel models. Second, an analysis describes how many data collectors per area 
are on average required to meet the outage probability for the given mean number of sensor 
nodes per area in case of Rayleigh fading. Third, this paper suggests a simple design method of 
the transmit power and the mean number of data collectors to meet the given outage probability. 
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The remainder of this paper is organized as follows: Section [TT] presents the system model 
based on a homogeneous Poisson point process (PPP). Section [III] analyzes the SIR distribution 
for Nakagammi-m fading channels and the SINR distribution for Rayleigh fading channels. 
Section [IV] derives the intensity of data collectors required to keep the outage probability below 
a certain value and suggests a design method of the transmit power. Section [V]discusses numerical 
results. Finally, Section [VI] concludes. 

II. System Model 

A sensor node senses or measures environments and then transmits its data to the closest 
data collector. Sensor nodes do not always have data to transmit but send them only when their 
sensing data are generated. For example, machines such as meters and event sensors may transmit 
data intermittently rather than continuously, and it has to be successful with probability above a 
certain value. In order to model intermittent transmissions, the sensor node's activity is defined 
as p. This value of p is between and 1, and this paper considers its small values. Meanwhile, 
data collectors that receive data from sensor nodes, are always ready to receive data from them. 

This paper considers environments where both of sensor nodes and data collectors are ran- 
domly deployed. Sensor nodes are distributed according to a homogeneous PPP, $ s , and they 
transmit sensed data to their nearest data collectors through random access schemes. \ s ,totai 
denotes the intensity of sensor nodes that is the average number of them per area. In order to 
consider unplanned deployments of data collectors, the random locations of data collectors are 
modeled as a homogeneous PPP, $ c , with intensity A c , like sensor nodes. Each sensor node 
transmits its data to a data collector closest to it, so a data collector builds a coverage based on 
Voronoi tessellation, as shown in Fig |2] 

The standard power loss propagation model with the path loss exponent a(> 2) and the 
Nakagami-m fading model are considered. In the Nakagami-m fading model |[22|. m = 1, 
m = (K + 1) 2 /(2K + 1) and m = oo model Rayleigh fading, Rician fading with parameter 
K, and no fading, respectively. Also, it is assumed that all sensor nodes transmit with the same 
power P. A typical data collector located on the origin receives the signal with Pr~ a Gs from 
a typical sensor node when the distance between them is r and the fading channel gain is 
Gs- By Slyvnyak's theorem [|23l . interfering nodes except for a typical sensor node located 
on X still constitute a homogeneous PPP with intensity X S) totai- Thus, the interference power 
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of the link between a typical sensor node and a typical data collector can be expressed as 
Ir — J2xj£<i> s \{Xa} P a Gij where Xj denotes the location of a interfering node and Gjj 
means the fading gain of a link between a typical data collector and a interfering sensor node 
j. {Gi,j} x e<b \{x } are independently and identically distributed (i.i.d.) random variables. Here, 
it is assumed that a typical data collector does not perform any scheduling for sensor nodes 
within coverage served by itself, so they may interfere with each other even though they are 
served by a common data collector. Eventually, the link of a transmitter-receiver pair experiences 
interference from interfering nodes distributed according to a homogeneous PPP with effective 
intensity A s = \ s> totai • p. If a sensor transmits using one of N resources that is chosen at random, 
\ s is X s ,totai • p/N. \ s decreases as N increases and this means that N is also a parameter for 
the system design. In this paper, N = 1 is assumed. 

When the interference is dealt with as noise and single antenna is equipped on both transmitters 
and receivers, the SINR is given by 

SINR = ^ 

\X \- a G s 



*s 



~ Ex, e * s \ { * 0} i^r"G /jj+ ^ 

where a 2 is the noise power and a 2 is equal to q p. In case of a 2 — > 0, © means the SIR. 

III. SINR Distribution 

This section analyzes the SINR distributions and derive simple-form SIR or SINR distribution 
for some specific channels. For more generalization of results, the fading gain distribution of © 
is first considered. 

Pr {G s > g} = ex P(-^) E a nk9 k (2) 

for some finite set J\f and a finite integer set /C. This type of complementary cumulative 
distribution function (CCDF) includes a variety of fading-gain distributions such as exponential 
distribution, chi-square distribution and gamma distribution. 

Lemma 3. 1 : Let sensor nodes and data collectors distributed with homogeneous PPP 's with 
intensities X s and A c , respectively and each sensor node builds communication link with a data 
collectors closest to it. When the CCDF of the fading gain of a desired signal is given by (|2|) 
and the fading gain of the interfering signal is denoted as a random variable, Gj, the CCDF 
of SINR is given by 
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Pr{SINR > 0} 

= 2vrA c E„eA^ Ekeic a n k (~(3) k (3) 

exp (— \ c nr 2 ) dr 



JO ' d( k 



(=nl3r a 

where a) = 7r(iT (l — E{Gf }, E{x} is the expectation of x andT{x) = / °° t x ~ x exp(—t)dt 
denotes the gamma function. The derivative in (TJ|) can be reexpressed as follow. 

d fc e X p(-A a 5(<»-C^) 



dC, k 

exp (-A s e(c, a) - c^ 2 ) Eto j E; =Q (-i) i+j G) w 

[A S C(C, a) + C^f & [A S C(C, a) + C^ 2 ] M 
Proof: See Appendix lAl ■ 
The result of SINR distribution in Lemma 13.11 requires cumbersome integrations and differ- 
entiations, but simple-form result can be obtained for specific channel models. 

To begin with, an analysis considers Nakagami-m fading channel. The received signal power 
experiencing Nakagami-m fading channel can be modeled using Gamma distributions. Thus, as- 
suming that desired signals experience Nakagami-m,; fading while interfering signals experience 
Nakagami-m, fading, the CCDF of fading gains can be give by 

Pr {G s > g} = £ L -f L e M-m s9 ) (5) 
fc=o K - 

mi_1 (ma) k 

Pr {Gj > g} = £ V-#L exp(-m^) (6) 
k=o 

© and © have the forms of (0, so the SINR distribution can be derived by using Lemma I3TI 
Generally, Lemma 13.11 requires the calculation of a derivative in © and it is too complex to 
calculate it for any m s , a and a. Fortunately, it is possible to obtain a simple form for the CCDF 
of SINR under interference limited environments, i.e. a 2 — > 0. 

Proposition 3.1: Let sensor nodes be randomly located with intensity \ s and served by the 
nearest data collectors randomly deployed with intensity intensity A c . When their links experience 
Nakagami-m fading given by © and ©), and a — > 0, the CCDF of SIR is given by 

Pr {SIR > 0} = , - Ac - 



k=0 k\ Ez=ol _i J ^k,l 



X c +\sC{mi,a)(m a p)a 

2 1 I \'J 

~T 

\c+XsC{mi,a)(m s f}) a 



where C(m, a) = m ~" T ^J {m+i) and A M = Ej =0 (-l) J ' ( •) U*=o ~ 1) ~ *] f° r 1 < k - 
Here, A ,o is defined as 1. 

Proof: See Appendix |Bj ■ 

In Nakagami-m fading model, m — 1 means the Rayleigh fading model. Thus, it is also easy 
to obtain the CCDF of the SIR for Rayleigh fading model. 

Corollary 3.1: Let sensor nodes be randomly located with intensity X s and served by the 
nearest data collector randomly deployed with intensity intensity A c . When all links experience 
Rayleigh fading with unit mean, and a — > 0, the CCDF of SIR is given by 



Pr {SIR > 0} = ^ T (8) 

\ c +X s C(l,a)ps 



where C(l,a) = T 1-i T (l + j- 



2 \ r< / 1 i 2 \ 2n 



a sin(27r/a) ' 

Proof: By substituting m s = 1 and m 8 = 1 into the results in Proposition |3J] ([8]) is obtained. 
Also, C(l,a) can be calculated by using the property of the gamma function T(l — z)T(z) = 



sin(7T2) ' 

As discussed before, when the noise power cannot be neglected, it is hard to obtain a simple 
form of the SINR for general m s because it requires the derivative of ©. However, when a path 
loss exponent, a, is equal to 4 and the fading channel is modeled as Rayleigh fading, the CCDF 
of SINR is simplified into a common integral form. 

Proposition 3.2: Let sensor nodes be randomly located with intensity X s and served by the 
nearest data collector randomly deployed with intensity intensity X c . when all links experience 
Rayleigh fading with unit mean and a path loss exponent a is 4, the CCDF of SINR is given by 

Pr{SINR > P} 

= 4^ exp ( ^\ K /: x A erfc ( (9) 

where K — ^- and erfc(x) = J^°exp (— t 2 ) dt is the complementary error function. 
Proof: The CCDF of fading gain for m s = 1 is the case of M — {1}, K, — {0} and a w = 1 
in ©• Thus, when a = 4 

Pr {SINR > p\ 

(10) 

= 27rA c / °°r ■ exp (-vrr 2 [A S /32C(1,4) + A c ] - /3r 4 a 2 ) dr 
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(flOl ) follows from © and (|23l ), and it can be evaluated by using the change of variables r 2 — > x 
and the integration formula, J* °° exp (— [ax + bx 2 }) dx = ^ exp fa^) erfc for a > and 

6 > 0. Here, K = vrC(l, 4) = ^. ■ 
Proposition 13.21 is the result for a = 4. When a is not 4, the CCDF of SINR can be expressed 
by generalized hypergeometric functions. But they are not simple, so this paper does not deal 
with them. 



IV. Intensity of Data Collectors 

When sensor nodes are spatially distributed according to a homogeneous PPP with a certain 
intensity, it is important to decide how many data collectors should be deployed in order to 
keep the success probability of random accesses above a certain value. This section analyzes the 
requirement of the intensity of data collectors deployed at random, and the effect of channels 
on its required intensity, given the intensity of sensor nodes and a target outage probability. The 
outage probability, e, is defined as Pr{SINR < (3 t } where f3 t is the minimal SINR value required 
for the successful receptions. 

The required intensity of data collectors for Rayleigh fading is presented in Corollary 14. II and 
Proposition 14.11 

Corollary 4.1: Let sensor nodes randomly located with intensity X s and served by the nearest 
data collectors. It is assumed that all links experience Rayleigh fading with unit mean and a — > 0. 
The necessary and sufficient condition of the intensity of data collectors randomly deployed, X c , 
for keeping the outage probability below e t , is 

\ c >±(l-e t )C(l,a)phs (ID 

where (7(1, a) is defined in Corollary I3.il 

Proof: (fTTI) can be directly derived from d8). ■ 
Proposition 4.1: Let sensor nodes randomly located with intensity \ s and served by the nearest 
data collectors. It is assumed that all links experience Rayleigh fading with unit mean and a = 4. 
The sufficient condition of the intensity of data collectors randomly deployed, X c , for keeping 
the outage probability below e t , is 



(1 - 2e t 



l + 8e t (l-e t ; 



(12) 



2 

where K = 

Proof: See Appendix O ■ 
The condition of A c in (fTTT) under interference-limited environments is necessary and sufficient 
while the condition in (fT2]) under environments with non-neglectable noise is just sufficient. In 
fact, (fT2l) has been derived from a lower bound of the complementary error function. But, for 
small values of a 2 , (fl~2)) also gives a tight lower bound of A c , and in particular, (fl~2T) is the same 
as (fTT) with a = 4 when a 2 — ► 0. 

Here, given (3 t , e t , a 2 and A s , the design method of the transmit power (P) of sensor nodes and 
the intensity (A c ) of data collectors is suggested, for a path loss exponent of four. The relations 
among these variables are given by ©, but it is not easy to use © directly for the design of P 
and A c . On the contrary, the lower bound of the CCDF of SINR with a simpler form of (|27T ) can 
give a simple design method for them. The lower bound of SINR CCDF in (T27T) is equivalent 
to the intensity condition of data collectors of (fl~2T ) and the second term within a square root 
in (fl~2)) approximately models the effect of noise. Now, the transmit power and the intensity 
of data collectors can be separately designed. First, for neglecting the noise effect, the second 
term within the square root of (fl~2T ) has to be much smaller that 1. The definition of a 2 gives a 
condition of the transmit power. 

Hp n f \ g2 / p < 8(j2 / p <r 1 (13) 

where (a) follows from the inequality of arithmetic and geometric means and the definition of 
K. Thus, the transmit power can be set to 

P = c% s (14) 

where c is a constant much less than one and it is a design parameter, c has to be set not 
only to neglect the noise power but also to keep transmit power as small as possible for sensor 
node's power saving. Next, the intensity of data collectors can be designed according to (fTTT ) 
because the intensity condition (fl~2l) is almost equal to (fTTT) if P is set by (fl~4l) . In this design, 
the transmit power is reciprocally proportional to \ 2 and this means that the longer the distance 
among sensor nodes is, the larger the required transmit power is, because of the noise effect, 
when the intensity of data collector is determined by ([TTb . Even though this design method is 
very simple, but it gives a good design method for the random deployment of data collectors to 
serve randomly distributed wireless sensors. Its performances will be shown in Section |V] 
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In interference-limited environments with Rayleigh fading channels, Corollary 14.11 shows the 
effect of the path loss exponents obviously. Because the function is a increasing function 
of < x < ir, it is obvious that the required density of data collectors decreases as the path loss 
exponent increases, when a > 2 and (3 t > 1 for given e t and A s , from the definition of (7(1, a) 
and CCD). 

On the contrary, it is not easy to express the required intensity of date collectors in case of the 
Nakagami-m fading with general m's, in a simple form. Here, the effect of wireless channels on 
system designs is analyzed by comparing the performances to those of Rayleigh fading, rather 
than deriving their requirements exactly, only when a 2 — > 0. When deploying data collectors 
with intensity A co for randomly distributed sensor nodes with intensity A s , let e(°) and 
denote the outage probabilities for the Rayleigh fading model and the another examined-fading 
model for the required SIR (3 t , respectively. First, in case of reference channel model assuming 
the Rayleigh fading, the A c o and e^°> have the following relation from ®. 

\ c , = 7 k(l-e^)C(l,a)pf\ s (15) 
On the other hand, in case of the examined-fading channel, A Cjm is defined as 

~X c , m = 7 k T (l-e^)C(l,a)phs (16) 

where is derived from ©. In other words, ([TBI means that the deployment of data collectors 
with A Cj0 in the Nakagami-m fading channel is equal to the deployment of data collectors with 
A Cjm in the Rayleigh fading channel in term of outage probability. Hence, A Cjm /A Ci0 quantifies 
the effect of wireless fading channels on the system design and is simplified from (fT5T ) and ([TBI ), 
as follows. 

\c,m _ £<°)(l-£ (m) ) _ l/e (m) ~l (-\H\ 

Xc,o ~ eW(l-eW) ~ 1/eW-l K J 

V. Numerical Results and Discussion 

This section evaluates and discusses the performance of systems with data collectors randomly 
deployed to serve randomly distributed wireless sensors, based on results of Section [Til] and 
Section [TV] It is assumed that the total intensity of sensor nodes (X s ,totai) spatially distributed 
according to a homogeneous PPP is lCT 2 m~ 2 . Also, p is set to 10~ 4 and it means that the sensor 
nodes awake on average every 1000 sec (about 17 minutes), when they transmit data to data 
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collectors during 100 msec on each awake mode. Also, the minimal SINR value ((3 t ) required 
for the successful reception of dB is considered. 

Fig. [3] shows the CDF of SINR according to P and y. This can be interpreted as the outage 
probability for (3 t which is a value on x-axis. P/a 2, s (or I /a 2 ) of 100 dB and 120 dB are 
assumed. These values mean that the transmit powers of sensor nodes are —10 dBm (0.1 mW) 
and 10 dBm (10 mW), when the power spectral density of the noise is —170 dBm/Hz and the 
bandwidth is 1 MHz. Fig. [3] indicates that analysis results in ([8]) and © definitely coincide 
with the simulation results. When P/a 2 is 100 dB, the ^'s of 10 and 20 result in the outage 
probabilities of 0.23 and 0.1, respectively. As P/a 2 increases, outage probability decreases. In 
other words, larger intensity of data collectors and higher transmit power lead to less outage 
probability. Fig. @] and Fig. |5] explain these effects more quantitatively. In Fig. H the outage 
probability decreases as the intensity of data collectors increases, and their required intensity 
can be obtained for a given outage probability. Also its lower bound by (|T2T) is shown. The lower 
bound of A c in (fT2l) is tighter when the effect of noise is reduced. The effect of noise on outage 
probability decreases as y- increases. This is because the increase of y- leads to the increase 
of received SNR because of the decrease in distances between data collectors and sensor nodes. 
Fig. [5] shows how the transmit power of sensor nodes affect the outage probability. The results 
of Fig. \5\ were evaluated by changing the intensity of sensor nodes for given relative intensities 
of data collectors. In other words, it shows the effect of noise by changing the geometric size 
of networks. The larger geometric size, i.e. larger distances between sensor nodes and data 
collectors, leads to the bigger effects of noise on the system performance. These results also 
verify that the design of the transmit power not only reduces the noise effect but also keeps 
the transmit power as small as possible. Also, under the environments of Fig. |3] the design 
by (TT4l) with c = 0.1 provides the transmit power of 9 dBm (i.e. P/a 2 = 119 dB), and it is 
observed that P/a 2 = 120 dB approaches the performance of P/a 2 — > oo in Fig.|3] These results 
confirm that (fl4l) is a very efficient design method. The path loss exponent is another crucial 
factor to have an effect on system performances. As Fig. [(^indicates, they result in very different 
performance for the same transmit power. At P/a 2 = 100 dB, the noise can be neglected in case 
of a pathloss exponent 3 while it causes severe performance degradation in case of a pathloss 
exponent 5. By contrast, when the noise effect can be neglected, larger a's result in less outage 
probability for given (3 t and A s . In fact, for given e t and A s , the required A c for a pathloss 
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exponent of a' increases by the factor of ^k^ Pt^ a » compared to a pathloss exponent of a 
when P/a 2 — > 0, where C(l, a) is defined in Corollary 14.11 For example, when f3 t — dB, the 
path loss exponents of 3 and 5 requires 1.54 and 0.84 times of the intensity of data collectors 
for the path loss exponent of 4. 

Fig. U}- Fig. [9] examine the performance for Nakagami-m fading channels. Fig. |7J explains how 
the line-of-sight factors of fading channels contribute to the SINR distribution. The increase in m 
results in the decrease in outage probability. But, m more than two does not have an big effect on 
the performance, compared to m equal to two. Fig. [7] also indicates that analysis results exactly 
coincide with simulation results when considering that the performance of P/a 2 = 120 dB is 
as good as that of P/a 2 — > oo. Fig. [8] shows the effect of channels on outage probability under 
the interference-limited environments. The outage probability decreases as m and the pathloss 
exponent increase. It means that the Rician fading and AWGN environments need less intensity of 
data collectors than the Rayleigh fading environments for the same path loss exponent. Moreover, 
from this figure, the intensity of data collectors required to meet a certain outage probability 
can be obtained. Fig. [9] examines the relative effect of other fading channels compared to the 
Rayleigh fading channel in term of the intensity of data collectors, which is defined in (TT7l) . 
it shows that m and a has a big effect on the system design such as the deployment of data 
collectors. 

So far, this paper analyzed and discussed the effect of the wireless channels, the transmit power 
and the intensity of data collectors on system performances when data collectors are randomly 
deployed to successfully collect the data from randomly-located sensor nodes. As the number 
of wireless nodes increases enormously in future, it is more and more difficult to design the 
system. For reducing these difficulties, efficient system design methods is required to deal with 
a huge number of wireless nodes, so the rigorous understanding about the spatial distribution 
and effect of interference will be basics for them. Even though this paper has considered only 
simple random access, these results will be able to be used as basic models for developing more 
sophisticated spatial resource management methods. 

VI. Conclusions 

This paper has considered the environment where receivers (data collectors) as well as trans- 
mitters (sensor nodes) are randomly deployed and each transmitter is served by the receiver 
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nearest to it. In network topology modeled by homogeneous Poisson point processes, analysis and 
simulation results showed the SINR distribution, and a simple design method of transmit power 
was suggested. Under interference-limited environments, the larger the path loss exponent and the 
portion of line-of-site factors were, the less the outage probability was. On the contrary, under 
non-neglectable noise environments, the large path loss exponent caused severe performance 
degradation. Moreover, the intensity of data collectors required to keep the outage probability 
above a certain value was derived, and it depends on required outage probability, an intensity 
of sensor nodes, a fading channel model, a path loss exponent and noise power. This required 
intensity helps to design such parameters as the amount of wireless resources and the access 
probability for medium access control. Random access scheme is very simple and does not cause 
control-overhead problems even under environments with a huge number of sensor nodes, but its 
required intensity of data collectors is never small. Thus, it is needed to find more sophisticated 
spatial resource management schemes and the result of this paper may be used as a basic model 
for them. 



This proof is similar to proof of theorem 1 in [19] that has considered the transmitter-centric 
coverage (or downlink) and only the transmitter intensity. Here, an analysis focuses on the 
receiver-centric coverage by data collectors (or uplink) and allows that multiple transmitters 
within the service area of a common data collector simultaneously transmit. For those differences 
and the completeness, this paper provides the full derivation of the CCDF of SINR. 

The probability that there is a data collector at a distance of r from a typical sensor node 
is 2tx X c dr. For this data collector to be a serving data collector of a typical sensor node, all 
other data collectors must be farther than r from a typical sensor node, and its probability is 
exp(— \ c 7cr 2 ). Thus, the probability density function of the distance between a typical sensor 
node and its serving data collector, f r (r), is equal to 2%\ c r ■ exp(— A c 7rr 2 ). 



Appendix A 



Proof of Lemma I3TT1 



The CCDF of SINR is 



Pr{SINR > (3} 




(18) 



27rA c / oo Pr{G' s > (3r a (I r + a 2 )} r exp(-A c 7rr 2 )dr 
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where I r = 'Ex j e^.\{x } \ x j\ * G i,j- From ©» 



Pr{G s > f3r a (I r + d 2 )} 

= Ej r {£n e Afexp(-n/3r Q [/ r + a 2 ])- 

E^an fc (/3r Q [J r + a 2 ]) fc } 
= EneM EfceK ank(fir a ) k - (19) 

E Ir {(I r + d 2 ) k exp(-n(3r a [I r + a 2 

(a) , ^ x fc ^E{exp(-C(J r +a 2 ))} 

— Z^neAf 2^k&K a nk\ — P' ) ^fc 



— l^neAf l^k£K a nk\—P r ) 



£=n/3r a 



where (a) and (b) follow from the definition of Laplace transform, £>x{() = Ex{exp(— (X)}, its 
•<x(t)(0 = (-!) fc ^ 



property, £ t k x(t)(() — (— l) fc d dc k ^ > an( * tne independence of I r and a 2 . The Laplace transform 



of L is 



£i r (0 = E /r {exp(-C/,)} 

= E 4>S)G/ {exp(-CEx i e* s \{^o} \ x j\ ° lG ij)} 
= E».{ILr i6 ».\{Jro}Eo«W(-CG/j I^D}} 
( = } exp (-2ttA s / °°[1 - E G/ {exp(-CG^- a )}]t;^) 



(20) 



(d) / o \ 

= exp (— 27rA s - 

/o°° (/o°°[l - exp(-(v~ a g}vdv) f Gl (g)dg) 
( = } ex P (-*^ir(-j)/ «^/ G/ (p)^) 

= exp(-A s £(C,a)) 

where (c) follows from the probability generating functional (PGFL) of the PPP ll23l : (d) uses 
the probability density function fdig) of a random variable Gf, (e) follows from the change of 
variable v~ a — > x and the definition of the Gamma function; (f) follows from the property of 
Gamma function xT(x) = T(l + x) and the definition of £(£, a). 
By substituting (OH) and <|2Q]) into (Ti8|) . © is derived. 

Also, © is obtained from the following equation which can be derived by the derivative of 
the exponential function and the chain rules: 



ft exp (/(*)) = 



exp(/W)Eto^E5=o(-l) J C)/0 



where r.) denotes .,,/ ! ... 



Appendix B 
Proof of Proposition |3JJ 

The fading gain of Nakagami-m fading channel given in ©) can be reexpressed as 

m s m a — l ^k 

Pr {G s >g}= exp(-n^) ]T — g k 



n=m a k=0 

So, the Nakagami-m fading is the case of J\f = {m s }, K. — {0, • • • ,m s 
Thus, 



1} and a nk 



£(£a) =7rCSr(l-| 

2 

m~ir( l-2-)r(m+2-' 



oo 2- g m i exp(— m^g) 
m i m T(mi) 



Jo 9 a 



dg 



where C(m, a) is defined as r 

When a 2 — > 0, the derivative © is calculated into 

d fc exp(-A s g(C,a)) 



exp (-\ s n(«C(m u a)) Eto J E^oM)^' (! 



j ^ r 



= exp (-A.ttC^CK, a)) Eto J E^oM)^' 

[-A s ^( fflll «)f [nto 1 (iii -j) - «)] c-' 

From (HH) and (ED, © is 



Pr{SIR > /?} 

(a) 



(b) 



2vrA c Er=o l! If(-/3) fc 

/ °°(m s /3)- fc Eto [AsTrCK, a)(m s /3)lr 
A fc)J • exp (-[A s 7rC(mi, a)(ra a/ fl)« + A c 7r]r 2 ) rdr 



m 8 -l (-l) fc 



27rA c E/£o 
v fc till 



fc! 



A s 7rC(mi,a)(m s /3) 



2l' 



A 



(c) 



/ °° exp (-[\ s 7iC(m t ,a)(m s (3)i + A c 7r]r 2 ) r 2 ' +1 dr 



m s -l (-1)* 



27rA c Efc=o 

^=0 /! 



A! 



2 

X s 7rC(mi, at)(m a p)« 



A 



k.i 



(d) 



-/-I 



A s 7rC(mi, a)(m s/ 9)« + A c 7r T(Z + 1) 



2" Z^i.=n n ' 
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where (a) follows from the definition of A k j in Proposition 13 . 1 1 (b) follows from the interchange 
of a summation and an integration, (c) follows from the calculation of the integral part by the 
definition of the Gamma function, and (d) follows from the property of the Gamma function 
T(l + 1) = l\ for a nonnegative integer I. 

Appendix C 
Proof of Proposition ED 



i 



L et T = K\c+Ki3 t a s an( j ^ __ — _ ^ e expressed as 



Pr{SINR > 0} = exp(r 2 )erfc(r)K 

2 r 
1+2t 



> exp(r 2 ) ■ 2 r exp( _ r2) . K (26) 



2 T K 



l+2r 2 

where (a) follows from the lower bound of the complementary error function. From © and (|26l . 

Pr{SINR>/3} >7rA " Ac+ ^ Aa 



2ftCT 2 +(7rA c +K/3 t 2 A s )2 (27) 

> I -e t 

where (b) follows from the definition of r and k. (1271) is rewritten into 

e t ^\l-[(l-2e t )vKphs]K 

(zo) 

-(i-^)(^ 2 AA 2 + 2A^ 2 ) >o 

which is a quadratic inequality with the form of a\ 2 c + 6A C + c > where a > and c < for 
< e t < 1. Thus, (1281) gives a positive lower bound of A c . By solving the inequality (T28T) for a 
variable A c > 0, (|T2]) is derived. 
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Fig. 2. Data collectors and sensor nodes distributed by a homogeneous Poisson point processes. Data collectors build the 
Voronoi tessellation (X S]t otai = 1CT 2 m~ 2 , A c = 5 x 10~ 3 m~ 2 , p = 0.01) 
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SINR (dB) 



Fig. 3. CDF of SINR for Rayleigh fading channels (a = 4; lines - simulation results; symbols - analysis results; star symbols 
represent the case of P/a 2 — > oo) 




Fig. 4. Effect of noise power on outage probability (a — 4; /3 t — dB; lines and symbols represent the exact performances 
and their approximations by (112) . respectively)) 
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Fig. 5. Effect of transmit power on outage probability (a — 4; fit = dB; solid circles and squares represent the transmit 
power values by the design of d!4t with c — 0.1 and c = 0.01, respectively) 




Fig. 6. CDF of SINR for Rayleigh fading channels according to the pathloss exponents (lines - simulation results; circles, 
squares and triangles represent the case of P/a 2 — > oo for a = 3, 4, 5, respectively) 
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Fig. 7. CDF of SINR for Nakagami-m fading channels (a — 4; lines - simulation results; symbols - analysis results when 

P/a 2 ->■ oo) 




Fig. 8. Effect of wireless channels on outage probability (fit = dB; P/a 2 oo) 
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